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Resumen 
En este trabajo se presenta la segunda versión del programa de 


cómputo DRENAS (drenaje agrícola subterráneo), el cual incluye tres 
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soluciones analíticas y una solución de elemento finito tipo Galerkin de 
la ecuación diferencial del drenaje agrícola en su forma unidimensional; 
puede ejecutarse bajo sistemas operativos Windows de 32 bits y 64 bits. 
Los modelos analíticos y numéricos incluidos en la versión 1.0 del 
programa se  reprogramaron en Visual Basic 2017, realizando 
modificaciones para mejorar sus capacidades de cálculo, registro de 
datos y visualización de resultados. El modelo numérico se amplió para 
resolver la ecuación de Boussinesq no lineal con coeficientes variables, 
incluyendo el término de recarga vertical, y para manejar las opciones 
de usar en la frontera de los drenes condiciones tipo Dirichlet, radiación 
lineal o radiación fractal. Los resultados proporcionados por el módulo 
de la solución numérica del programa fueron validados 
satisfactoriamente haciendo uso de una solución analítica y un 


experimento de drenaje. 


Palabras clave: ecuación de  Boussinesq, coeficiente de 
almacenamiento variable, condición de radiación fractal, propiedades 


hidráulicas del suelo. 


Abstract 


The second version of the computer program DRENAS (subsurface 
agricultural drainage) is presented, which includes four analytical 
solutions and a finite element solution Galerkin type of the differential 
equation of agricultural drainage in its one-dimensional form and can be 
run under Windows systems of 32-bit or 64-bit. The analytical and 
numerical models included in version 1.0 of the program were migrated 
from Visual Basic 2005 to Visual Basic 2017 and modified to expand and 
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improve their calculation, data recording and results visualization 
capabilities. The numerical model can now solve the non-linear 
Boussinesq equation with variable coefficients including the vertical 
recharge term and allows use at the boundary of the drains, Dirichlet 
condition, linear radiation or fractal radiation. The results provided by 
the module of the numerical solution of the program were validated 


using an analytical solution and a drainage test. 


Keywords: Boussinesq equation, variable storage coefficient, fractal 


radiation condition, soil hydraulic properties. 


Recibido: 14/03/2018 
Aceptado: 29/07/2019 


Introducción 


Los terrenos agrícolas que presentan problemas de salinidad y/o de 
manto freático somero pueden ser recuperados o controlados mediante 
drenaje artificial que en el caso de la parcela puede ser del tipo 
subterráneo. La variables a determinar en el diseño de un sistema de 
drenaje subterráneo es la profundidad de los drenes, la separación de 
los mismos y su diámetro. El cálculo de estas variables debe considerar 
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las características del suelo (textura, estructura, propiedades físicas e 


hidráulicas) y la condición del flujo de agua en el subsuelo. 


Es posible realizar estudios de los procesos de transferencia de 
agua en los sistemas de drenaje agrícola con la ecuación de Boussinesq 
de los acuíferos libres, que aunque considera de forma simplificada las 
transferencias de masa y energía en la zona vadosa del suelo, 
proporciona descripciones generales del flujo de agua en el espesor 
saturado del medio poroso (Ritzema, 2006). Para realizar modelaciones 
detalladas de la dinámica del agua con la ecuación de Boussinesq, en el 
caso de acuíferos libres someros se debe considerar la dependencia del 
coeficiente de almacenamiento del acuífero respecto de la carga 
hidráulica y usar la condición de frontera en los drenes que mejor 
represente el flujo de drenaje. Por una parte, Fuentes, Zavala 4 
Saucedo (2009) establecen formalmente la relación entre la curva de 
retención de humedad y el coeficiente de almacenamiento en acuíferos 
libres y, por otra, Zavala, Fuentes 8 Saucedo (2007) demuestran que la 
transferencia de agua del suelo hacia el interior de los drenes debe ser 


descrita con una condición de radiación no lineal. 


Zavala, Saucedo € Fuentes (2014) desarrollan la primera versión 
computacional para aplicar la ecuación de Boussinesq con coeficiente de 
almacenamiento variable sujeta en los drenes a condiciones de radiación 
fractal, a la cual denominan DRENAS (drenaje agrícola subterráneo). Sin 
embargo, la herramienta de cómputo se programó en Visual Basic 2005, 
por lo que opera sólo en sistemas Windows de 32 bits, y sus opciones 


gráficas y base de datos dependen de Microsoft Office 2003. 
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El objetivo general de este trabajo fue desarrollar en Visual Basic 
2017, la segunda versión del programa de cómputo DRENAS, para 
expandir sus capacidades de cálculo, mejorar su interfaz gráfica, y 
eliminar su dependencia de programas externos para el manejo de 
datos y elaboración de gráficos. En la actualización también se 
incorporaron y programaron nuevas opciones de cálculo, a fin de 
incrementar las alternativas de simulación, como el manejo del término 
de recarga o descarga de la ecuación de Boussinesq como una función 
dependiente del tiempo; representación de la condición inicial de la 
carga hidráulica como una función variable en el espacio; posibilidad de 
usar restricciones mecanicistas para los parámetros de forma de la 
curva de retención de humedad de Van Genuchten (1980) cuando se 
utiliza para definir el coeficiente de almacenamiento como función de la 
carga hidráulica, y la opción de usar condiciones de frontera en el dren 
del tipo Dirichlet. Con la meta de depurar el programa de cómputo, se 
realizó su validación, considerando una solución analítica para drenaje 
subterráneo en régimen estacionario y datos de una prueba de drenaje 


experimental de laboratorio. 


Materiales y métodos 


Ecuaciones de base 
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En el programa de cómputo se resuelve la forma completa o simplificada 
de la ecuación de Boussinesq unidimensional del drenaje agrícola. En 


general, se puede escribir como: 
UM E= 2 [TME| + Ryo (1) 


Donde H = H(x, t) es la carga hidráulica contada a partir de un 
estrato impermeable o nivel de referencia (L) y es función de la 
coordenada horizontal x y del tiempo t; T(H) es la transmisibilidad del 
acuífero (L?T*); para el caso de acuíferos libres es T(H) = K.H; K: es la 
conductividad hidráulica a saturación del suelo (LT*); u(H) es el 
coeficiente de almacenamiento del acuífero; Ry(t) es el volumen de 
recarga o descarga en la unidad de tiempo por unidad de área del 


acuífero (LT?). 


La recarga o descarga Ry(t) es una variable difícil de determinar, 
ya que depende de las condiciones de flujo existentes en la superficie de 
suelo (lluvia, evaporación, tirante de agua), el flujo del agua y el 
régimen de humedad en la zona vadosa, la estratigrafía de la zona no 
saturada del suelo, la macroporosidad, grietas, etcétera. En esta nueva 
versión del modelo de simulación numérico se incorporó el término 


Rw(t), representándolo mediante un polinomio que depende del tiempo: 
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Ry (t > 0) = a, t? + b,t? +c,t + d, (2) 


Donde a,, b,, Cr y d, son coeficientes a calcular, considerando 
mediciones experimentales o estimaciones que se tengan de la 


evolución en el tiempo de la recarga para un evento específico. 


Se consideró en el programa la relación entre el coeficiente de 
almacenamiento y la curva de retención de humedad establecida por 
Fuentes et al. (2009): 


u(H) = 65 — 01 — Hs) (3) 


Donde 8; es el contenido volumétrico de agua a saturación (L?L”?y; 
y Hs es la elevación de la superficie del suelo (L). Es posible obtener una 
representación analítica explícita del coeficiente de almacenamiento a 
partir de la Ecuación (3) si la curva de retención de humedad del suelo 
es conocida. Se seleccionó la relación clásica de Van Genuchten (1980), 
ampliamente aceptada en estudios de campo y de laboratorio por su 


flexibilidad descriptiva: 


8(p) = 0, + (05 — 0,111 + (Y/Ya)"]” (4) 


Donde 8, es el contenido volumétrico residual (L3L”%); ya < 0 es un 
parámetro de escala del potencial de presión del agua en el suelo v, (L); 
m y n son parámetros de forma adimensionales. En esta nueva versión 
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del programa de cómputo se incluyeron relaciones mecanicistas entre m 
y n que satisfacen la teoría de la infiltración, la restricción clásica de 
Burdine (1953) m = 1 - 2/n, así como las relaciones fractales de 
Fuentes, Brambila, Vauclin, Parlange 8 Haverkamp (2001), que derivan 
en su estudio de la conductividad hidráulica de los suelos no saturados. 
Las restricciones fractales son las denominadas poro neutro m = (1 — 
4s/n)/s, poro geométrico m = (1 - 2s/n)/s y poro grande m = (1 - 
4s/n)/2s; donde s es la dimensión cociente del objeto fractal, tal como 
se considera el suelo en ese estudio, siendo s la razón entre la 
dimensión fractal del objeto D; (Mandelbrot, 1982; Falconer, 2014) y la 


dimensión del espacio de Euclides E, s = Df /E. 


Al introducir la Ecuación (4) en la Ecuación (3) se obtuvo la 


siguiente relación: 


u(H) = (05 0,41 — [1 + (4H —= Hs) /Pay13 (5) 


Se retuvo en esta segunda versión del programa la condición de 


radiación fractal para los drenes derivada por Zavala et al. (2007): 


qa = YKinlUa — Do) / PY? (6) 


Donde qa es el flujo de drenaje; y es un coeficiente de 
conductancia adimensional; Kin = vWKsKa) es la conductividad de la 
interfaz suelo-dren; Ks es la conductividad hidráulica a saturación del 
suelo; Ky la conductividad de la pared del dren; P, la profundidad de los 
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drenes; Ha, la carga hidráulica en la posición del dren; Do, el espesor del 
acuífero; s = (1/2)(s¡+S4) es la dimensión cociente en la interfaz suelo- 
dren; s, es la dimensión cociente del suelo y sy es la dimensión cociente 
de la pared del dren. La relación entre la porosidad areal (pea) y la 


porosidad volumétrica (p) se obtiene de acuerdo con la idea probabilista 


COMO area, = $ “p*= p*, se tiene que la porosidad volumétrica: 


(1-p) +6 =1 (7) 


y la porosidad areal: 


1 
1 = 
(1 — Marcar): + sn =1 (8) 


La conductividad de la pared del dren se obtiene con una fórmula 
basada en la ley de Poiseuille: Ky = (1/2)(9/W)Harea, a(Ruo)?, donde g es 
la aceleración de la gravedad; v, la viscosidad cinemática del agua; Ryp, 
el radio hidráulico del dren. Se introdujo en el versión 2.0 de DRENAS la 
opción de manejar directamente la condición de radiación lineal (5 = 0.5 
en la Ecuación (6)) sin necesidad de utilizar las rutinas de cálculo de s; y 


Sa: 


Modelaciones numéricas simplificadas del drenaje agrícola se 
pueden llevar a cabo imponiendo en los drenes condiciones de frontera 
tipo Dirichlet (carga hidráulica en la posición espacial del dren). En la 
nueva versión del programa de cómputo se incorporó esta condición de 
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frontera como una alternativa a la condición de radiación (6); se 
seleccionó la siguiente función, que depende de la variable 


independiente tiempo: 


rat > 0) = Uat + bat? + Cg + dut (9) 


Donde Xdren es la coordenada horizontal donde se ubica el dren; y 
da, Da, Ca Y da son coeficientes que se deben calcular a partir de 
mediciones que se tengan de la evolución de la carga hidráulica sobre el 
dren. La Ecuación (9) tiene la flexibilidad de describir comportamientos 


extremos de la carga sobre el dren, como el abatimiento y ascenso. 


Por último, para realizar la simulación numérica con la Ecuación 
(1) debe definirse el estado inicial de la carga hidráulica en el sistema. 
En este trabajo se incluyó la opción de usar hasta un polinomio de tercer 


grado, que es función de la coordenada horizontal x: 


Hí(xt= 0) ==, +bixé +0, +0 (10) 
p p p p 


Donde ap, Dp, Cp Y dp son coeficientes a determinar a partir de 
mediciones de la carga hidráulica en el sistema en el instante inicial de 


cálculo o simulación. 


En la solución numérica del sistema (1)-(10) se empleó el método 
del elemento finito tipo Galerkin para la discretización espacial; un 


esquema de diferencias finitas para la discretización temporal; el 
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método iterativo de Picard para la linealización del sistema resultante, y 
un método de gradiente conjugado precondicionado de para la solución 
del sistema de ecuaciones algebraicas (Zavala et al., 2014; Zienkiewicz, 
Taylor, 4 Zhu, 2013; Noor 8 Peters, 1987). El sistema de ecuaciones 


resultante de la discretización se programó en Visual Basic 2017. 


Es conocido que cálculos simplificados de separación entre drenes 
y abatimiento del manto freático pueden ser realizados aplicando 
soluciones analíticas obtenidas para formas reducidas de la ecuación de 
Boussinesq (Ecuación (1)); por ejemplo, para régimen de flujo 
permanente se tiene la fórmula de Hooghoudt (1940), y en régimen 
transitorio las relaciones de Glover-Dumm (Dumm, 1954) y Fuentes et 
al. (1997), las cuales ya se incluían en la versión original del programa 


de cómputo y se retienen en la actualización. 


Desarrollo de la interfaz gráfica 


En la segunda versión del programa de cómputo DRENAS, se 
reprogramaron en Visual Basic 2017 todos los módulos de captura de 
información de la versión original; se desarrollaron bases de datos 
internas independientes de Microsoft Access; se desvinculó de Microsoft 
Excel el tratamiento de las gráficas, y se generó el archivo ejecutable 
para instalar DRENAS 2.0 en cualquier computadora que disponga de 
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sistema operativo Windows de 64 bits, e incluso se generó también la 


opción para sistemas Windows 32 bits. 


El programa DRENAS 2.0 tiene dos módulos de cálculo: uno 
denominado soluciones analíticas y otro modelo numérico (Figura 1). El 
módulo soluciones analíticas tiene cuatro secciones: dos para realizar 
cálculos para condiciones de flujo de agua en régimen permanente 
usando la solución de Hooghoudt (1940) y dos para condiciones de flujo 
en régimen transitorio; uno para realizar cálculos con la solución de 
Dumm (1954) y el otro con la solución de Fuentes et al. (1997). En 
estas cuatro secciones se puede calcular la separación entre drenes o el 
módulo de drenaje. Respecto a la versión original del modelo DRENAS, 
en la versión 2.0 estas secciones fueron mejoradas, al incorporar las 
siguientes opciones: a) almacenamiento de las simulaciones realizadas 
en una base de datos interna; b) carga de ejemplos previos; Cc) 
generación de un reporte de la simulación y alternativa para exportarlo 
a un archivo pdf; d) La gráfica de resultados también puede ser 
exportada como imagen. Un ejemplo de una sección de cálculo del 


módulo de soluciones analíticas se presenta en la Figura 2. 
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Figura 1. Pantalla general del programa DRENAS 2.0. 
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Caiga hdrautcas (him 


os 


25 


Era dl a tn. ts: De 


Figura 2. Una sección de cálculo del módulo soluciones analíticas. 


El módulo de la solución numérica se programó en una forma Tab 
Index (en inglés), que contiene seis pestañas o secciones: las primeras 
cuatro están desarrolladas para la captura de datos referente al sistema 
de drenaje, características del dren, propiedades del suelo y datos 
necesarios para la simulación numérica, como son los relacionados con 
el control de tiempo, discretización espacial, condiciones límite y 
visualización gráfica de resultados; la quinta pestaña sirve para 
almacenar los datos del ejemplo en una base de datos interna, y la 


última pestaña contiene el botón para ejecutar la simulación numérica. 
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Las nuevas capacidades del módulo de simulación numérica de la 
versión 2.0 de DRENAS que se programaron se establecieron de la 


siguiente forma: 


a) Primero de la base de datos de suelos UNSODA 2.0 (Nemes, 
Schaap, Leij € Wósten, 2001) se seleccionaron las curvas 
experimentales de retención de humedad de 208 suelos reportados en la 
literatura; se ajustaron con la Ecuación (4), considerando las 
restricciones mecanicistas entre m y n antes descritas (Zavala, Saucedo, 
8: Fuentes, 2018), y los resultados obtenidos se incluyeron en una base 
de datos interna del programa de cómputo para ponerlos a la 
disposición de los usuarios del mismo (Figura 3). Si se selecciona la 
opción de usar datos de esta base del programa, al elegir un suelo, las 
secciones de propiedades físicas, hidráulicas y coeficiente de 
almacenamiento se llenan automáticamente con sus parámetros 
correspondientes. La mejora a resaltar en esta nueva versión es que se 
pueden manejar las relaciones entre m y n de los modelos de Burdine, 


poro neutro, poro geométrico y poro grande (Figura 4). 
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TAR Y 


Figura 3. Módulo numérico: base UNSODA ampliada. 
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Modelo para la curva de mención de ver, Garucitas (19001 
Corterado volunálnco a anturación (E) 


Corterido volumétnco ressdual (Or) 


Pecbrmetro de escaín (9d) 
Parámetro de forma (n) 
Perámatro de farma (m) 


E Ps ! = e E 1 E 


Figura 4. Módulo numérico: propiedades suelo-agua. 


b) En la solución numérica de la ecuación de Boussinesq se 
incorporó el término de recarga Ry(t) como una función del tiempo 
(Ecuación (2)); se realizaron las modificaciones de programación 
necesarias para poder imponer en los drenes la condición de frontera 
tipo Dirichlet (Ecuación (9)), y se programó también la condición inicial 
tipo polinomio cúbico (Ecuación (10)). En la Figura 5 se presenta la 


nueva sección desarrollada para incluir estas tres condiciones. 
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Drenaje Agrícola MO 
Unidades y Sistema de drenaje E Salir 


»4 Modelo Numérico 


Sist. de Drenaje (datos completos) | Dren (Datos completos) | Prop. del Suelo | Simulación | Guardar Ejemplo | Ejecutar 


Discretización | Registro de resultados | Condiciones limte 


R(0) = a, + b,t? + c,t + dr 
(0<x<L) 


087 
(A 


Figura 5. Módulo numérico: condiciones límite y recarga. 


c) Se desarrolló una nueva sección para registrar los datos de la 
simulación, cargar datos para comparar simulaciones y para seleccionar 
el tipo de gráfica a visualizar en tiempo real que no depende de 


programas externos (Figura 6). 
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Figura 6. Módulo numérico: gráficas de la simulación. 


Resultados y discusión 
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Los resultados del programa de cómputo se compararon contra 
resultados externos (solución analítica y una prueba de drenaje en 
laboratorio) para detectar y corregir errores de programación, así como 


revisar la consistencia de las soluciones codificadas. 


Validación 1 


Se seleccionó la solución analítica del tipo Hooghoudt presentada en 
Fragoza et al. (2003), que considera radiación lineal (s = 0.5 en la 
Ecuación (6)) en los drenes para describir la evolución de la superficie 
libre en un sistema de drenaje subterráneo con recarga constante (Ro) y 


flujo de agua en régimen permanente: 


H?(0) =H¿+(Ro/KIx(L —x) (11) 


donde H. es la carga hidráulica al centro entre drenes. Si se define 


h(x,t) = H(x<,t) - Do y considera la condición de radiación lineal se tiene: 


ha = [y (+ y)?D2 +8(Q + y)h¿(h¿ + 2D,) — (4 + y)Do]/[2(2 + y)] (12) 
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donde hy es la carga hidráulica sobre el dren (x=0 y x=L); D, es el 
espesor del acuífero; y y el coeficiente de conductancia de la condición 
de radiación lineal que se usa en los drenes. Analizando la Ecuación (2) 


se tiene para este caso que a,=b,=C,=0 y d,=Ro. 


Los datos del sistema de drenaje son los reportados por Fragoza 
et al. (2003), derivados de un experimento de campo realizado en el 
distrito de riego 076 Valle del Carrizo, Sinaloa. El sistema de drenaje 
evaluado tiene una separación entre drenes L = 50 m, profundidad 
media de los drenes P = 1.5 m, profundidad del estrato impermeable Do, 
= 3.5 m, elevaciones del estrato impermeable y de la superficie del 
suelo H;=0.0 m y Hs=5.0 m. La parcela tiene suelo de textura arcillosa, 
con una conductividad hidráulica a saturación de K; = 0.557 m/d, y el 
coeficiente de almacenamiento u = 0.1087. Se asume y = 1.5 y h¿= 0.5 
m. Con la Ecuación (12) se tiene ha = 0.365 m y el valor de recarga 
obtenido con la condición para régimen permanente donde la descarga 
por el dren es igual a la recarga, Ro = 0.000944 m/d. Con estos valores 
se aplicó la Ecuación (11) para obtener la distribución espacial de la 


carga hidráulica a lo largo de la separación entre drenes. 


Se reprodujo con el programa de cómputo el escenario descrito, 
modelándose el proceso transitorio flujo de agua en el sistema hasta 
alcanzar el régimen estacionario, el cual corresponde al descrito con la 
solución analítica. Se asumió la condición inicial de carga hidráulica 
constante H(x, t=0) = 4.50 m (por tanto, ap = Dp= Cp= O y d, = 4.50 m 
en la Ecuación (10)) y para tener las mismas condiciones de la solución 
analítica, en el programa se usó la condición de frontera de radiación 


lineal, 5=0.5 en la Ecuación (6). Se consideró la hipótesis de igualdad 
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entre la conductividad del suelo y la de la pared del dren K; = Kg, lo cual 
genera Kin = Ks = Kg. Se reinterpretó el valor del coeficiente de 
conductancia reportado por Fragoza et al. (2003), para adaptarlo a la 


forma de radiación que se maneja en DRENAS 2.0, obteniéndose y = 


0.045. La discretización espacial se realizó con el número de nodos N 
1000 (Ax = 0.05 m), y en la discretización del tiempo se usó Ati; = 
1.157E-05 d, Atmax = 6.94E-04 d, Y Atmin = 1.157E-06 d. El tiempo de 


simulación fue de 720 días. 


En la Figura 7 (a y b) se presentan los resultados de la evolución 
del manto freático a la semana y al año de drenaje libre. Se observa que 
la solución numérica tiende a la solución analítica conforme avanza el 
tiempo y la reproduce de forma correcta. Esta comparación es una 
evidencia de que los módulos de captura de información de la solución 
numérica del programa y en particular el módulo de cálculo que maneja 
el término de la recarga vertical, condición inicial y condición de 
frontera, operan bien y están libres de errores de programación; 
también se observa que las series simuladas son estables y libres de 


oscilaciones numéricas. 
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Tiempo (dias) Carga al centro entre drenes (m) Tiempo (dias) Carga al centro entre drenes (m) 


7 18784975219815 4 AT ADA TA 65 287849751939 4 0003244303904 


a) Primera semana de drenaje b) Un año de drenaje 


Figura 7. Simulación de un escenario de drenaje agrícola (coeficiente 


de almacenamiento constante y recarga vertical constante). 


Validación 2 


Zavala, Fuentes 8 Saucedo (2004) realizaron un experimento en un 
módulo de drenaje en el que se instalaron dos drenes de 30 cm de 
longitud (/p) y 5 cm de diámetro (Dp), con N, = 233 perforaciones 
circulares de diámetro d, = 0.158 cm, distribuidas de modo uniforme en 
la superficie de cada tubo. Las dimensiones de cada dren: L = 100 cm, P 
= 120 cm y Do = 25 cm. El módulo se rellenó con suelo de textura 


arenosa, que se saturó aplicando una carga de agua constante en la 
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superficie del suelo hasta eliminar el aire. Con los drenes tapados se 
removió el exceso de agua en la superficie y se cubrió para eliminar la 
evaporación. Por último, se destaparon los drenes y se midió durante 10 
días (240 horas) el volumen de agua drenado; se debe notar que la 
condición inicial corresponde a H(x,0) = P + Do (dp, = 145 cm en la 
Ecuación (10) y el resto de sus coeficientes son nulos) y la recarga es 


nula (Ry = O en la Ecuación (2)). 


La porosidad del suelo fue O, = 0.5396 cm*/cm”, y con este valor 
se resolvió la Ecuación (7), obteniéndose si; = 0.7026. El valor de la 
conductividad hidráulica a saturación reportado es K; = 18.3 cm/h. El 
coeficiente de almacenamiento del acuífero descrito con la Ecuación (5) 
para los parámetros 6; = Oy, 0, = 0, Wa = -41.8 cm y n = 2/(1-m) = 
3.19 (restricción de Burdine, 1953). 


La porosidad areal del dren eS Harea, y = Ao/Aa = 0.0097 cm*/cm? y 
al resolver la Relación (8) se tiene sy = 0.5688. Con el radio hidráulico 
de los orificios circulares de la pared del dren Ryp = 0.0395 cm, la 
viscosidad cinemática del agua v = 36 cm“/h se aplica la ley de 
Poiseuille para calcular la conductividad de la pared del dren Ky¿ = 
2,670.8 cm/h. Los valores de la conductividad en la interfaz suelo-dren 
y la dimensión relativa son Kin = 221.08 cm/h y s=0.6357. El 
coeficiente y de la radiación fractal reportado por Zavala et al. (2004) se 
reinterpretó para adaptarlo a la forma de radiación que se maneja en 
DRENAS 2.0 (Ecuación (6)), obteniéndose y = 0.0749. 


El coeficiente de almacenamiento variable y la condición inicial de 
saturación total del medio poroso originan que la simulación deba ser 
ejecutada empleando pasos de tiempo pequeños al inicio de la misma, 
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lo cual minimiza problemas de convergencia y estabilidad numérica. Los 
pasos de tiempo utilizados son At»; = 2.77E-04 h, Atin; = 2.77E-05 h y 


Atmax = 0.2 h. El número de nodos para discretizar el espacio es de 200. 


En la Figura 8 (a y b) se presenta, respectivamente, el proceso de 
captura de datos de la simulación y la evolución de la lámina drenada. 
En la Figura 8b se puede observar que el programa de cómputo 
reproduce de manera satisfactoria los datos experimentales de lámina 
drenada, teniéndose una solución numérica estable y libre de 
oscilaciones, lo cual es un indicador de que el módulo de coeficiente de 
almacenamiento variable trabaja de forma correcta. Los pasos de 
tiempo y espacio utilizados en la simulación son adecuados, ya que la 
lámina drenada total calculada (23.915 cm) difiere de la lámina medida 
(23.92 cm) en sólo 0.03%. 


«7 EmmGrafica DES 
3 
en 
0 48 96 a 192 240 
Pa 
a] Tiempo (hr) Lámina (cm) | Guardar gráfica_| 
240 DNISOSIS5 TES | Generar reponso_| 
< > Preso de IA O 
a) Captura de información: b) Láminas drenada 
coeficiente de almacenamiento variable experimental y simulada 


Figura 8. Simulación de un escenario de drenaje agrícola (coeficiente 
de almacenamiento variable, recarga nula y condición de radiación 
fractal en los drenes). 
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Conclusiones 


La nueva versión del programa de cómputo DRENAS (versión 2.0) es útil 
para describir el funcionamiento hidráulico de drenes subterráneos, para 
diseñar un sistema de drenaje y resolver problemas inversos de 
caracterización hidrodinámica de suelos a partir de eventos de drenaje. 
Esta versión de la herramienta computacional incorpora las siguientes 
funciones: a) bases internas para almacenar los datos ingresados en 
cada uno de los ejemplos que se analizan; b) base de suelos UNSODA, 
ampliada con parámetros para los modelos de conductividad hidráulica y 
retención de humedad del poro neutro, poro de la media geométrica y 
poro grande; c) opción para representar la recarga o descarga del 
acuífero como una función del tiempo; d) opción para representar la 
condición inicial de la carga hidráulica como función del espacio y usar 
condiciones tipo Dirichlet en los drenes variables en el tiempo; e) genera 
reportes de las simulaciones y permite exportar las gráficas de 
resultados. La visualización gráfica de las modelaciones evita la 
dependencia de programas externos, y se ejecuta sobre sistemas 


operativos Windows de 32 y 64 bits. 
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